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ABSTRACT 

We present scaling laws for advection, radiation, magnetic drag and ohmic 
dissipation in the atmospheres of hot giant exoplanets. In the limit of weak 
thermal ionization, ohmic dissipation increases with the planetary equilibrium 
temperature (T eq > 1000 K) faster than the insolation power does, eventually 
reaching values > 1% of the insolation power, which may be sufficient to in- 
flate the radii of hot Jupiters. At higher T eq values still, magnetic drag rapidly 
brakes the atmospheric winds, which reduces the associated ohmic dissipation 
power. For example, for a planetary field strength B — 10 G, the fiducial scaling 
laws indicate that ohmic dissipation exceeds 1% of the insolation power over the 
equilibrium temperature range T eq ~ 1300-2000 K, with a peak contribution at 
T eq ~ 1600 K. Evidence for magnetically dragged winds at the planetary ther- 
mal photosphere could emerge in the form of reduced longitudinal offsets for the 
dayside infrared hotspot. This suggests the possibility of an anticorrelation be- 
tween the amount of hotspot offset and the degree of radius inflation, linking the 
atmospheric and interior properties of hot giant exoplanets in an observationally 
testable way. While providing a useful framework to explore the magnetic sce- 
nario, the scaling laws also reveal strong parameter dependencies, in particular 
with respect to the unknown planetary magnetic field strength. 



1. Introduction 



Hot giant exoplanets, including hot Jupiters, are among the best characterized exoplan- 
ets. Secondary eclipses, transmission spectroscopy and orbital phase curves in particular 
have provided constraints on the atmospheric properties of several transiting members of 
this class (see Charbonneau 2009; Winn 2010 for reviews). Increasingly detailed radiative 
and circulation models for the atmospheres of these exoplanets have also been developed to 
help interpret the rich emerging phenomenology (see Burrows & Orton 2010; Showman et 
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al. 2010 for reviews). Finally, evolutionary models have been systematically used to infer 
the bulk interior properties of hot giant exoplanets with precisely measured transit radii (see 
Baraffe et al. 2010 for a review). 

It was recently proposed that magnetic effects could have important consequences for 
the atmospheric and bulk interior properties of hot giant exoplanets. Magnetic drag and 
ohmic dissipation, which result from kinematic induction by thermally-driven winds in the 
weakly-ionized atmospheres of hot giant exoplanetsjl] will indeed brake atmospheric winds 
(Perna et al. 2010a; Rauscher & Menou 2011) and deposit extra heat in the deep atmosphere 
or the interior adiabat (Batygin & Stevenson 2010; Perna et al. 2010b), thus affecting some 
of the observable properties of these exoplanets. 

Here, we establish simple scaling laws for the strength of magnetic drag and ohmic 
dissipation in the atmospheres of hot giant exoplanets, as a function of their radiative equi- 
librium temperature. These scaling laws allow us to study the nature of the transition into 
the magnetized regime, as the radiative equilibrium temperature is raised above ~ 1000 K, 
and to explore how the limit of strong coupling with the planetary magnetic is eventually 
reached, at higher temperatures. The remainder of this paper is organized as follows. In 
§2, we recall the basic mechanism of magnetic induction in the atmosphere of a hot giant 
exoplanet. In §3, we establish scaling laws for the strength of magnetic drag and ohmic dis- 
sipation as a function of the planetary radiative equilibrium temperature. The main results 
emerging from these scaling laws are described in §4. In §5, we discuss various parameter 
dependencies, the strength of induced magnetic fields, connections to previous work on this 
topic and a few plausible observational signatures of magnetic effects in the population of 
hot giant exoplanets. We conclude in §6. 



2. Basic Mechanism 

Day-side insolation on a hot Jupiter generates fast atmospheric winds (see, e.g., reviews 
by Showman et al. 2008; 2010). As the weakly-ionized gas flows across the planetary 
magnetic field, which is presumably anchored to the planet's bulk rotation via deep-seated 
electric currents in the convective interior, an additional magnetic field is induced in the 
atmosphere. In steady-state, this magnetic induction is balanced by resistive diffusion and 
associated ohmic dissipation in the resistive atmosphere. 



1 We are excluding from our analysis giant exoplanets which are hot because of their young age. While 
magnetic effects may also be important for them, our focus is on older planets with atmospheric temperatures 
and winds determined by strong external irradiation, as exemplified by hot Jupiters. 
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Let us consider a spherical coordinate system (r, 9, <p) that is rotating with the planet. 
Let us further assume that the planetary magnetic field takes the form of an aligned dipole 
(of surface strength, B^ p ) and that the dominant atmospheric flow is zonal (azimuthal) in 
nature: 3> V r , Vg. To leading order, the induced current J satisfies the steady-state 
resistive induction equation 



where rj is the local resistivity of the weakly-ionized atmospheric gas. Note the implicit 
assumption that the deep-seated electric currents generating the dipolar planetary field are 
located well below the atmospheric region of interest, so that they can be omitted from the 
above induction equation (Liu et al. 2008; Batygin & Stevenson 2010; Perna et al. 2010a,b). 

For a purely zonal flow and a strictly aligned dipole, a purely toroidal field {B^) is 
induced, which is sustained by purely poloidal currents: J = (c/47r)V x (see, e.g., Liu et 
al. 2008). Magnetic induction does not operate at the rotational equator for an aligned dipole 
because vertical speeds are negligibly small and meridional wind speeds do not contribute. 
The aligned dipole assumption is likely to remain qualitatively correct even for moderate 
levels of misalignment. 

Focusing on the dominant, zonal component of the atmospheric flow, a representative 
horizontal momentum balance for the steady flow in the azimuthal (e^) direction may be 
written as 



where the zonal acceleration (LHS) resulting from the day-night pressure gradient (first term 
on the RHS) is reduced by the bulk Lorentz force due to magnetic drag (second term on the 
RHS; e.g., Zhu et al. 2005). 

Our main goal with this work is to study how this momentum balance changes with the 
atmospheric resistivity, which depends itself strongly on the planet's radiative equilibrium 
temperature via thermal ionization balance. Considering the above induction and momen- 
tum equations together, one anticipates the following behavior. At low enough atmospheric 
temperatures (high resistivity rj), magnetic drag (oc J) will not significantly reduce the zonal 
wind speed, (in Eq. [2j). At fixed and B dip values, the induced current J will thus 
increase in inverse proportion to the resistivity rj (in Eq. PQ). At large enough temperatures 
(low resistivity rj), however, the winds will be strongly dragged and it is no longer obvious 
how wind speeds and induced currents should scale with resistivity in Eq. ([1]). It is thus of 
great interest to establish the range of atmospheric temperatures over which such a transition 
occurs and to clarify what happens in the limit of high conductivity (high temperatures). 




(1) 




-4- 



3. Magnetic Scaling Laws 

Our goal in establishing magnetic scaling laws for hot giant exoplanet atmospheres is 
not to obtain an accurate description of magnetic effects for this class of planets, as this will 
require detailed mult i- dimensional models. Rather, our focus is on the leading-order scaling 
with the radiative equilibrium temperature, which sets the atmospheric resistivity. As we 
establish these scaling laws, we will need to invoke important simplifying assumptions. 



3.1. General Properties 

We adopt fiducial parameters for the scaling laws that represent a typical hot Jupiter, 
with a hydrogen-dominated, solar composition atmosphere, so that a perfect gas constant 
1Z = 4.59 x 10 7 erg/g/K and a specific heat at constant pressure C p = 1.43 x 10 8 erg/g/K are 
appropriate. The planetary radius is R p = 10 10 cm, with a surface gravity g = 890 cm s~ 2 . 
We explore surface magnetic field values from £>di P = 3 G to £>di P = 30 G. Planetary radia- 
tive equilibrium temperatures ranging from T eq = 1000 K to 2200 K are considered, which 
correspond to semi-major axes a ~ 0.09-0.02 AU around a Sun-like star Jf| 



3.2. Model Atmospheres 

For the purpose of the present study, we focus our discussion on the most actively forced, 
"weather layer" region of the atmosphere, defined here as the region located between the 
thermal (infrared) photosphere above and the insolation (visible) photosphere below, which 
corresponds to a level above which ~ 50% of the insolation flux has been absorbed (see details 
below). For simplicity, the vertical thermal structure of the atmosphere is approximated with 
the one-dimensional radiative solution discussed by Guillot (2010; see also Hansen 2008) for a 
zero incidence angle and a dilution factor / = 0.5 representing a dayside average. As shown 
by Guillot (2010; see his Fig. 4), this solution provides reasonable temperature-pressure 
profiles for hot Jupiter atmospheres over a range of orbital radii, even though it clearly is a 
simplified approach and it neglects important three-dimensional aspects of the atmospheric 
structure. In what follows, we refer to this pressure- dependent temperature profile for the 
dayside weather layer as Td ay - Note that, by construction, Td ay > T eq everywhere in this 



1/2 

2 The radiative equilibrium temperature is denned by T cq = T» for zero Bond albedo, where 

and i?* are the stellar effective temperature and radius, and D is the planet-star orbital separation (e.g., 
Guillot 2010). 
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weather layer model. 

Following Guillot (2010), a single thermal absorption coefficient = 10~ 2 cm 2 g _1 
and a single visible absorption coefficient k s = 4x 10~ 3 cm 2 g -1 are adopted. For these 
values, the thermal and visible photospheres (defined by r = 2/3) are located, respectively, 
at P t h ~ 60 mbar and P v ~ 150 mbar (above which ~ 50% of the stellar flux has been 
absorbed). The number of pressure scale heights over which the weather layer extends 
vertically in our models is thus modest, AlnP = hiP v — lnP t h — 0.9. An internal heat flux 
corresponding to T int = 150 K is adopted for the thermal structure calculation, but the exact 
value has little effect on the atmospheric properties of the weather layer of interest here. 

Guillot (2010) has also shown that these values for the grey absorption coefficients are 
adequate to reproduce the global radiative equilibrium balance of the prototypical hot Jupiter 
HD209458b, with T eq ~ 1420 K. The fixed values of the absorption coefficients adopted in 
our scaling laws, despite the large range of T eq = 1000-2200 K considered, is arguably one of 
the strongest simplifying assumptions we make in our models. We comment on the effects 
of variations in opacities and other model parameters in §5.11 

3.3. Ionization Balance 

Thermal ionization, rather than photoionization, is expected to be dominant over the 
range of atmospheric pressures (P > 50 mbar) of interest for the scaling laws (see, e.g., Perna 
et al. 2010a; Batygin & Stevenson 2010). Rather than using the same approximate solution 
to the Saha equation as in Perna et al. (2010a,b), which is only justified at low enough 
temperatures (T < 1800 K), we use here a more general but still simplified solution to the 
Saha equation (e.g., Sato 1991). In the weakly-ionized regime, when only single ionization 
has to be considered, the ionization fraction x e (<^C 1) satisfies the equations 

Tig r > Tlj . . 

j 

where n e and n n are the number densities of electrons and neutrals, respectively (in 
cm -3 ), rij is the number density of element j, Ij is the corresponding first ionization potential, 
T is the temperature, n is the total number density of the gas and the other symbols have 
their usual meaning. Solar composition for the first 28 elements of the periodic table (H to 
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Ni) is assumed for the ionization calculations. In the limit x 3 ; <C 1, which ceases to be valid 
for some elements at T > 2000 K for the densities of interest here, the above expression 
essentially reduces to that used by Laughlin et al. (2011). Based on this solution for the 
ionization balance, the local electric resistivity of the atmospheric gas is evaluated everywhere 
in the weather layer as rj = 230VT/x e cm 2 s" 1 (see Perna et al. 2010a for details). 

The ionization solutions obtained with Eq. (j4j) are still approximate, however, since all 
the degeneracy factors in the Saha equation have been arbitrarily set to unity. This crude 
simplification is justified by the fact that the specific atmospheric composition of any given 
hot Jupiter is a priori unknown and could easily deviate from solar in a non-trivial way. In 
addition, from the point of view of the scaling laws, inaccuracies in the degeneracy factors 
are dwarfed by the exponential dependence with temperature appearing in Eq. (jJJ), given 
the wide range of temperatures considered. Indeed, we verified that factor several changes in 
the degeneracy factors entering the Saha equation have little practical impact on our scaling 
law results. 



3.4. Forcing, Momentum Balance and Dragged Wind Speeds 

Showman et al. (2010) describe how to obtain a simple, order-of-magnitude scaling for 
the zonal wind speed in a hot Jupiter atmosphere, in the absence of magnetic drag, based 
on steady non-linear balance between the zonal acceleration term and the pressure gradient 
term in the horizontal momentum equation. The same dimensional analysis of Eq. fl2]), 
taking the limit B^p = 0, leads to the same scaling 

Vj, = ^Tl AT horiz AlnP, (4) 

where ATh or i z is the typical day-night temperature differential along the equator, 1Z is 
the perfect gas constant and A In P measures the number of vertical pressure scale heights 
in the weather layer over which this horizontal temperature differential applies. 

In the presence of magnetic drag, a similar dimensional analysis of Eq. (JSJ leads to the 

VI _ TZAT hoTiz AlnP V^B 2 

R p R p Airpq' [ ' 

where the bulk Lorentz force has been evaluated via dimensional analysis of the steady-state 
resistive induction equation (Eq. [T|) for the latitudinal current: Jq ~ cV ^B / \4nrrj) (see, e.g., 
Perna et al. 2010a). We obtain an order-of-magnitude estimate of the zonal wind speed in 
the presence of magnetic drag by solving for the positive root of this quadratic equation for 
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An additional complication with the above estimate is that the day-night temperature 
differential itself, AT horiz , depends on the zonal wind speed via the efficiency of heat 
advection from the dayside to the nightside. This results in a non-linear coupling between 
advection and radiation processes in the atmospheric flow. While Showman & Guillot (2002) 
proposed to capture this dependence with an exponential relation between AThoriz and the 
ratio of the atmospheric advective and radiative timescales, we adopt a power-law dependence 
here for simplicity, keeping the index as a free parameter. The horizontal temperature 
differential between the day and the night sides of a tidally-locked, hot giant exoplanet is 
thus estimated as 

AT horiz = min ^ (^) 1 , (6) 

A I V T rad/ 

so that when r a d v /Trad > 1, the full differential Td ay /2 is applied, while if r a d v /r ra d < 1, 
advection reduces the effective temperature differential between day and night in proportion 
to the values of these two timescales. A default value n = 1 is adopted for the power law 
index but other values are explored below. 

The advective timescale is evaluated as 

r adv = tt, (7) 
'<t> 

while the radiative timescale is evaluated as 

where g is the gravitational acceleration, P and Td ay are the pressure and dayside tem- 
perature at the location of interest in the weather layer and a is the Stefan-Boltzmann 
constant (e.g., Goody & Yung 1989). 

Equations (jSJ) and (jSJ) for our two unknowns, and AThoriz, are solved iteratively. 
As it turns out, for the majority of the parameter space of interest here, atmospheres are 
strongly radiative so that the precise form of the scaling with T a d v /T ra d or n in Eq. fl6]) has 
little consequences for our main conclusions. 

The ohmic power dissipated per unit volume in relation to the magnetic drag is evaluated 

as 



4 71T]J 2 _ Vp?_ 
C 2 Alt 7] 



Qohm ^ a ' (9) 
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Integrating vertically, and then horizontally over the volume of the weather layer, this 
yields an estimate of the total ohmic power dissipated as a result of induction in the resistive 
atmosphere, 

P ohm = 4vrW Q ohm x H p d\nP, (10) 

J AlnP 

where the local pressure scale height is evaluated as H p = TZT day /g. It is worth empha- 
sizing that this estimate of the ohmic power is not expected to be accurate at better than 
an order of magnitude level, given that it overlooks important issues related to the geometry 
and the spatial distribution of electrical currents within and outside the weather layer (see, 
e.g., Liu et al. 2008; Batygin & Stevenson 2010; Perna et al. 2010b). 



4. Results 

Figure [1] shows our estimate for the zonal wind speed, at the thermal photosphere 
(Pth — 60 mbar) as a function of the planetary radiative equilibrium temperature, T eq . 
Results for an assumed magnetic field strength B d i P = 3, 10 and 30 G are shown as solid, 
dashed and dash-dotted lines, respectively. 

These curves cover three different regimes for momentum balance in the thermally- 
forced atmospheres of hot giant exoplanets. At low T eq values, advection reduces the day- 
night temperature differential ATh or i z according to Eq. (JSJ), which leads to the relatively 
steep slope of with T eq . As shown best by the solid line, at higher T eq values, this slope 
becomes shallower as the atmosphere becomes radiatively dominated (when r adv /r rad > 1), 
with a negligible role for advection (as defined by Eq. [6]). Finally, the three curves show 
that magnetic drag eventually becomes efficient at braking the winds and that it does so 
at lower radiative equilibrium temperatures for stronger magnetic field strengths. Results 
comparable to those shown in Fig. [T] are obtained for deeper levels in the weather layer than 
the thermal photosphere. 

The exponential decline in V$ at high T eq values, which is caused by strong magnetic 
drag, can be understood as follows. As resistivity drops, the dominant balance in the mo- 
mentum equation (Eq. [3]) eventually becomes one where the pressure gradient acceleration 
term is balanced by the deceleration term from the bulk Lorentz force. In that limit, 
becomes small and simple dimensional analysis in that limit shows that it scales in propor- 
tion to the resistivity, oc rj, which itself is an exponentially decreasing function of the 
atmospheric temperatures in the weather layer, and thus of T eq . 
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Fig. 1. — Estimate of the zonal wind speed, (in km/s), at the thermal photosphere as a 
function of the planetary radiative equilibrium temperature, T eq (in K). Solid, dashed and 
dash-dotted lines show results for surface magnetic field strengths -Bdip = 3, 10 and 30 G, 
respectively. The exponential decline of at high T cq is caused by strong magnetic drag. 
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Fig. 2. — Ratio of advective to radiative timescales, r adv /r rad , at the thermal photosphere 
as a function of the planetary radiative equilibrium temperature, T eq (in K). Solid, dashed 
and dash-dotted lines show results for surface magnetic field strengths -Bdip = 3, 10 and 
30 G, respectively. The radiative equilibrium temperature above which advection becomes 
negligible (r a d v /Trad > 1) is comparable to the temperature at which the strong magnetic 
drag regime begins for B dip ~ 30 G. 
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Fig. 3. — Ohmic power (in Watts) dissipated as a result of induction in the weather layer as 
a function of the planetary radiative equilibrium temperature, T eq (in K). Solid, dashed and 
dash-dotted lines show results for surface magnetic field strengths -Bdip = 3, 10 and 30 G, 
respectively. Ohmic power rises exponentially with T eq until the regime of strong magnetic 
drag is reached, at which point the power starts declining exponentially. 
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Fig. 4. — Ratio of ohmic to irradiation power (P hm/-Pirr) as a function of the planetary 
radiative equilibrium temperature, T eq (in K). Solid, dashed and dash-dotted curves show 
results for surface magnetic field strengths -Bdip = 3, 10 and 30 G, respectively. Assuming that 
a ratio -P hm/-Pirr > 1% is required for sizable planetary radius inflation, this figure illustrates 
how radius inflation should only exist over a restricted range of equilibrium temperatures, 
according to our magnetic scaling laws. 
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Figure [2] illustrates the nature of the change from an advective to a radiative behavior 
for hot giant exoplanet atmospheres as the radiative equilibrium temperature, T eq , is raised 
above ~ 1200 K. Below 1200 K, r adv /r rad < 1 and advection reduces the day-night temper- 
ature differential according to Eq. (jH]). Above 1200 K, however, the atmospheric thermal 
structure (at the thermal photosphere) becomes essentially radiative, with little role, if any, 
for advection. As shown by the three curves in Fig. [2j this transition is effectively acceler- 
ated by magnetic drag, when it significantly reduces the wind speeds, the more so for strong 
magnetic field strengths. 

Figure [3] shows the ohmic power P hm dissipated as a result of induction in the weather 
layer as a function of the planetary radiative equilibrium temperature, T eq . The exponential 
rise and subsequent exponential decline of P hm with T eq can be understood as follows. In 
the limit of weak magnetic drag, is independent of the resistivity and it varies only weakly 
with T eq . Equation (jUJ) then implies that P h m oc Qohm oc 1/rj. On the other hand, in the 
limit of strong magnetic drag, since cx T], the ohmic power declines exponentially with T eq 
according to P h m oc Qohm oc V- As shown by Fig. [3j the radiative equilibrium temperature 
at which the ohmic power peaks depends sensitively on the value of the surface magnetic 
field strength. 

Figure H] shows a different representation of the ohmic dissipation results, where the 
ohmic power P hm has been scaled to the total irradiation power P irr = 47r Pp<rT e 4 q received 
by a hot giant exoplanet with radiative equilibrium temperature, T eq . The ratio Pohm/Prr is 
useful because it has been suggested that a value > 1% is needed for sizable radius inflation 
of a hot Jupiterjfl by deposition either in the bulk convective interior (Batygin & Stevenson 
2010) or in the deep atmosphere (Perna et al. 2010b). Our simple magnetic scaling laws lack 
the descriptive power to identify where in the planet the ohmic power is dissipated, but they 
clarify the range of radiative equilibrium temperatures over which such power is in principle 
available for radius inflation. Like the peak ohmic value, the range of radiative equilibrium 
temperatures over which Pohm/Prr > 1% (or any other fixed ratio) is a sensitive function of 
the surface magnetic field strength. 



3 We note that any energetic threshold for radius inflation, such as the ratio -P hm/-firr, will generally be 
a function of the planet's mass and age (e.g., Miller ct al. 2009), so that the 1% value adopted here should 
be taken as merely indicative. 
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5. Discussion 

5.1. Parameter Dependencies 

An examination of the various assumptions made in deriving our magnetic scaling laws 
makes it clear that they can only capture the most basic, order-of-magnitude behavior for 
hot giant exoplanet atmospheres across the wide range of radiative equilibrium tempera- 
tures considered. It is important for the validity of the scaling laws that the central physical 
ingredient of the theory, atmospheric resistivity, varies exponentially with atmospheric tem- 
peratures, via thermal ionization balance. Still, there is some freedom in renormalizing the 
scaling laws by adopting different parameter values than the fiducial ones we have adopted. 
It is thus important to understand the sensitivity of the scaling law results to variations in 
the different model parameters. 

We find that modest changes (~ 20%) in the assumed values of the planetary radius, 
R p , the gravitational acceleration, g, the perfect gas constant, 71, or the specific heat at 
constant pressure, C p , have only minor quantitative effects on the scaling law results. Raising 
the assumed gas metallicity to ten times solar, which increases the free electron density 
accordingly, reduces the temperature at which the strong drag regime is entered (from, 
e.g., ~ 1600 K to ~ 1475 K, for B^ p = 10 G), but the scaling law results are otherwise 
qualitatively similar to the solar metallicity case. Similarly, changing the power law index 
used to evaluate the day-night temperature differential in Eq. ([6]) from n = 1 to n = 3 or 
n = 5 further steepens the slope of V$ with T eq at low T eq values (leftmost region in Fig. [Q, 
but the scaling law results remain otherwise qualitatively similar. 

The scaling law results are sensitive to variations in model parameters which directly 
affect the atmospheric temperatures in the weather layer, however, as expected from the 
exponential dependence of the resistivity with temperature. For example, increasing the 
thermal opacity coefficient, fc t h, by 50% nearly doubles the ohmic power and shifts the ohmic 
peak from ~ 1600 K to ~ 1500 K, for = 10 G. Reducing k t ^ by a factor of two leads to 
about three times less ohmic dissipation, a peak ohmic power at ~ 1650 K for B^ p = 10 G 
and magnetic Reynolds numbers reduced by a factor of a few. Increasing the visible opacity 
coefficient, k v , by 50% reduces the ohmic power by a factor ~ 5, shifts the ohmic peak to 
slightly lower temperatures and reduces magnetic Reynolds numbers by a factor of a few. 
Reducing the value of k v by a factor of two increases the ohmic power by a factor ~ 5, shifts 
the ohmic peak to slightly larger temperatures and increases magnetic Reynolds numbers by 
a factor of a few. 

Similarly, an increase in the assumed dayside dilution factor, /, entering Guillot's ra- 
diative solution raises the ohmic power somewhat and shifts its peak to lower temperature 
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(from, e.g, ~ 1600 K for / = 0.5 to ~ 1475 K for / = 0.7, for B dip = 10 G). A reduction in 
the magnitude of the maximum day- night temperature differential in Eq. ([6]) from 0.5 T day 
to 0.3 T day reduces the ohmic power by about a factor two and the peak values by ~ 25%. 
Conversely, an increase of this maximum day-night temperature differential from 0.5 T<} ay to 
0.7 Tday doubles the ohmic power and increases the peak value by ~ 20%. 

Another factor which may significantly impact atmospheric temperatures in the weather 
layer, and thus the magnetic scaling laws, is the possibility of a temperature inversion in 
the upper atmosphere of strongly-irradiated giant exoplanets (e.g., Fortney et al. 2008). 
Hot giant exoplanets with such inversions would typically have reduced temperatures in 
their weather layer, from extra absorption of stellar light at altitude. Even though we did 
not explore this scenario in any detail, and instead adopted a simple set of Guillot (2010) 
radiative solutions with > n v for our thermal profiles, one can anticipate a systematic 
shift of the magnetic scaling law results at fixed T eq for those planets with inversions, in 
proportion to the temperature deficit in their weather layer. A detailed investigation of 
the effects of temperature inversions on magnetic scaling laws would constitute a valuable 
extension of our work. 

Besides the strong dependence on the planetary magnetic field strength, which is illus- 
trated explicitly in Figs. 1-4, our parameter space exploration reveals a strong dependence of 
the scaling law results on parameters affecting the atmospheric temperatures in the modeled 
weather layer. This suggests that improved treatments of the atmospheric thermal structure, 
e.g. based on three-dimensional modeling, would greatly benefit efforts to better understand 
the magnetic behavior of hot giant exoplanet atmospheres as a function of their radiative 
equilibrium temperature. 



5.2. Magnetic Reynolds Number and Induced Fields 

As the radiative equilibrium temperature of a hot giant exoplanet is raised, the degree 
of coupling between the atmospheric flow and the planetary magnetic field also increases. 
The nature of this coupling in the strong drag limit identified by our scaling laws, at the 
high end of the range of equilibrium temperatures considered, is an important issue. 

Figure shows an estimate of the magnetic Reynolds number of the atmospheric flow, 
R m , as a function of the planetary radiative equilibrium temperature, T eq , using the same 
notation as in previous figures. This magnetic Reynolds number is evaluated as 

Rm = (11) 
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where is the zonal wind velocity at the thermal photosphere (Fig. [I]), H p is the pres- 
sure scale height and r\ is the resistivity. Fig. [5] reveals that the magnetic Reynolds number 
is less than unity and independent of the magnetic field strength at low enough radiative 
equilibrium temperatures. Beyond T eq ~ 1300 K, however, R m exceeds unity and reaches 
different, near-asymptotic values for different strengths of the planetary magnetic field. The 
near-asymptotic behavior of R m at high temperatures can easily be understood as result- 
ing from the previously established scaling oc rj in the strong drag limit, with different 
normalizations for the different magnetic field strengths. 

Traditionally, the regime R m > 1 would indicate a significant degree of coupling be- 
tween the flow and the magnetic field, and the possibility of sustained dynamo action in 
the atmosphere. It is unclear, however, whether the moderately large values of R m found 
in Fig. [T] are sufficient to trigger a self-sustained dynamo, and what the nature of such a 
dynamo might be, given the unusual setup of a very shallow atmospheric flow threaded by 
the planetary magnetic field. 

Irrespective of the issue of sustained dynamo action, the significant gas-field coupling 
arising from values of R m > 1 raises the possibility of large induced field strengths even 
in the simple kinematic framework adopted in this study. In particular, a strong enough 
induced toroidal magnetic field (B^) could in principle result in secondary induction of a 
strong poloidal field component by meridional atmospheric motions. This could challenge 
one of our basic assumption, which is that a dominant contribution to the poloidal field 
component is provided by the planetary dipole field (-Bdip m Eq. [T]). 

Some insight on this issue may be gained by simple dimensional analysis. Using the 
same scaling as in §3.41 for the latitudinal current, Jg ~ cV ^B^j '(Anrf), and assuming that 
the relevant gradient scale is ~ H p in the relation between the induced current and the 
induced field, Jg = (c/47r)V x B^, one deduces a simple estimate of the strength of the 
induced toroidal field in terms of the background dipole field: B^ ~ R m B<n p . In the strong 
drag limit, this suggests values for the induced toroidal field which can be in excess of the 
planetary dipole field. 

If large enough field values were induced into a poloidal component by meridional mo- 
tions, Eq. ([1]), which is a starting point for our scaling laws, could cease to be valid: the value 
of the poloidal field may no longer be set by the fixed value of the planetary dipole field but 
rather by multi-dimensional induction in the atmospheric flow itself. It is unclear, however, 
whether such large induction of poloidal field is expected in the weather layer of hot giant 
exoplanets. Indeed, simple dimensional analysis analogous to the one performed above sug- 
gests an induced poloidal field strength B po \ ~ R^B^, where R* m = VgH p /r] is the magnetic 
Reynolds number associated with the meridional (rather than zonal) flow. Therefore, to the 



-17- 




Fig. 5. — Logarithm of the magnetic Reynolds number, R m , of the atmospheric flow at the 
thermal photosphere, as a function of the planetary radiative equilibrium temperature, T eq 
(in K). Solid, dashed and dash-dotted lines show results for a surface magnetic field strength 
-Bdip = 3, 10 and 30 G, respectively. The exponential drop in caused by strong magnetic 
drag is responsible for the near-asymptotic values of R m reached at high T eq . 
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extent that the meridional flow velocities (Vg) are much weaker than the zonal ones, in the 
presence of separate drags acting on each component, may still be small enough for the 
induced poloidal field to remain negligible compared to the background dipole field. 

Our scaling law results in the strong drag limit assume that the planetary dipole field 
remains dominant relative to any induced poloidal field. It is difficult to evaluate if or when 
this assumption might break down without considerably more detailed explorations of the 
physics of magnetized atmospheric flows than performed here, such as mult i- dimensional 
models with distinct drag treatments for the zonal and meridional flows, or perhaps even 
full MHD treatments to explore dynamo issues. Until such studies become available, it 
will thus be important to remain cautious when using simple scaling arguments like ours to 
interpret the phenomenology of hot giant exoplanet atmospheres in the strong drag limit. 

5.3. Connection to Previous Work 

Our scaling law results are broadly consistent with previously published results for 
the zonal wind speeds in the weather layer of hot giant exoplanets. In particular, the 
values shown in Fig. [1] for zonal wind speeds are qualitatively consistent with the zonally- 
averaged results of atmospheric circulation models for the hot giant exoplanet HD209458b 
(T eq ~ 1400 K) with weak (£>di P = 3 G), moderate (B dip = 10 G) and strong (-B<j ip = 30 G) 
imposed drag, as reported in Perna et al. (2010a). Extrapolating the trend at low T eq values 
in Fig. [1] also points to zonal wind speeds in broad agreement with those reported by Lewis 
et al. (2010) for the atmosphere of the hot Neptune GJ436b (T oq ~ 650 K) and it supports 
the notion that the atmospheric flow on this cooler exoplanet can effectively be treated as 
free of magnetic drag. 

The nature of the transition from the unmagnetized to the magnetized atmosphere 
regime for hot giant exoplanets is a particularly interesting issue to examine in light of our 
scaling law results. In terms of magnetic drag, Fig. [T] indicates that the magnetized regime 
starts around T eq ~ 1200-1600 K, depending on the strength of the planetary magnetic field. 
Defining the value of T eq above which ohmic dissipation begins to be important is more 
difficult because it depends on a threshold value for the ratio of ohmic to irradiation power 
above which significant radius inflation is expected. This ratio is not well defined because it 
depends on details of the ohmic dissipation process but a threshold value of -P hm/-Pirr ~ 1% 
is probably a reasonable guess (Guillot & Showman 2002; Batygin & Stevenson 2010; Perna 
et al. 2010b). As shown in Fig. |U this 1% threshold is reached at T eq ~ 1100-1500 K 
according to our scaling laws, with a strong dependence on the planetary magnetic field 
strength. When evaluating the transition in the magnetized regime, one should remember 
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the great degree of simplification used in deriving our scaling laws, which is exemplified by 
some of the strong parameter dependencies discussed in §5.11 The nature of the spatial 
distribution of ohmic dissipation will also likely influence the transition temperature for the 
onset of radius inflation (Batygin & Stevenson 2010; Perna et al. 2010b). 

In their study of the heavy element content of giant exoplanets, Miller & Fortney (2011) 
suggest a temperature of ~ 1000 K as the threshold below which radius inflation appears to 
be absent from their sample of exoplanets. This temperature threshold is broadly consistent 
with, but somewhat lower than, the threshold values inferred on the basis of our scaling 
law results, shown in Fig. HJ It may be possible to bring these threshold values into closer 
agreement by lowering the critical value of P hm/-Pirr required for sizable radius inflation 
(below our 1% assumption) and/or by postulating that planetary magnetic fields in the 
exoplanet sample considered by Miller & Fortney (2011) are generally strong (I?dip ^30 G). 

In an analysis sharing some similarities with the present one, Laughlin et al. (2011) 
have studied a correlation between a measure of radius inflation for hot giant exoplanets, the 
so-called radius anomaly, and their effective temperature^ in an observationally constrained 
sample. The emergence of a positive radius anomaly at T e g > 1000 K, as shown in their 
Fig. 2, is broadly consistent with our scaling law results. It is unclear whether or not our work 
supports the power law dependence of the radius anomaly with T e s advocated by Laughlin 
et al. (2011). This relation depends on the intricate physics that connects ohmic dissipation 
with the effective radius inflation and it simply is not addressed by our scaling laws. It is 
conceivable, however, that the large spread in the amount of radius anomaly found at any T e g 
value (see their Fig. 2) could be accounted for by variations in the strengths of the planetary 
magnetic field in the sample, since we find -Bdip to be a rather sensitive parameter of the 
theory. We note that there also is an interesting indication of a systematic drop in radius 
anomaly at T e g > 1900 K in the Fig. 2 of Laughlin et al. (2011). It is possible that this 
feature is associated with the drop in ohmic dissipation, and thus in the amount of radius 
inflation, that is expected in the strong drag regime according to our scaling law results (see 
our Fig. @J. 

Batygin et al. (2011) have studied in detail some of the consequences of ohmic dissi- 
pation for the thermal evolution of hot giant exoplanets. In their analysis, these authors 
assume that the global efficiency, e, of ohmic heating relative to that from stellar irradiation 
is a model constant, fixed at values of 1, 3 or 5%. Our scaling law results challenge this 
simplifying assumption since we find that the ratio -P hm/-Pirr is a strong function of both 



For strongly-irradiated planets with negligible internal heat fluxes, effective and radiative equilibrium 
temperatures are equivalent. 
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the planetary effective temperature, T e g, and the planetary magnetic field strength (see our 
Fig. HJ). The emergence of sizable radius inflation at T cS > 1400 K in the evolutionary 
models of Batygin et al. (2011; see their Figs 6-8) is certainly consistent with the typical 
threshold values suggested by our scaling laws, but the existence of a well-defined transition 
temperature would require a rather uniform set of planetary magnetic strengths according 
to our analysis. Given the difference in methodology, we note that it may be possible to 
recast some of the results of Batygin et al. (2011) in the language of our scaling laws, by 
considering that their assumed ohmic efficiency values of e = 1, 3 and 5% correspond to 
different effective planetary magnetic field strengths. 

Batygin et al. (2011) argued that magnetic drag has a very limited effect on the magni- 
tude of zonal winds in hot giant exoplanet atmospheres. This is inconsistent with our scaling 
law results, which do suggest the existence of a strong drag regime, as shown in our Fig. (TJ 
and it also conflicts with the strongly dragged simulation results described in Perna et al. 
(2010b) and Rauscher & Menou (2011). Our work suggests that the emergence of a strong 
drag regime is essential to the saturation and the eventual decline of ohmic dissipation at 
high planetary effective temperatures. It is worth noting that a systematic reduction in the 
efficiency of ohmic dissipation e in the models of Batygin et al. (2011) when T c g > 1600- 
1900 K would likely bring their model radii in closer agreement with the data (see their 
Figs. 6-8). Such a drop of e with T e g could also mitigate the emergence of unstable, over- 
flowing evolutionary tracks identified by Batygin et al. (2011) for models with high enough 
values of T cS (> 1400-1800 K). 



5.4. Observational Signatures 

Magnetic drag and ohmic dissipation are separate manifestations of the same magnetic 
induction process operating in the atmosphere of a hot giant exoplanet. It is of great interest 
to explore the type of observational signatures expected from magnetic drag and ohmic 
dissipation since these observables should ultimately be related to each other. As discussed 
earlier, the most direct observable manifestation of ohmic dissipation is planetary radius 
inflation. Signatures of magnetic drag, on the other hand, would most likely emerge in the 
form of reduced atmospheric wind speeds. 

Over the past few years, it has become possible to constrain the magnitude of wind 
advection in the atmosphere of a few hot giant exoplanets, by measuring their infrared 
phase curves (e.g., Harrington et al. 2006; Cowan et al. 2007; Knutson et al. 2007, 2009a,b; 
Crossfield et al. 2010). The standard interpretation that has been given to an infrared peak 
emission that is offset from the planet's orbital phase is that heat is advected away from the 
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dayside by eastward equatorial winds at the planet's thermal photosphere (e.g., Showman 
et al. 2008, 2009, 2010; Rauscher & Menou 2010; Dobbs-Dixon et al. 2010; Burrows et al. 
2010; Heng et al. 2011; Showman & Polvani 2011). Magnetic drag will act to reduce the 
speed of zonal winds and thus the amount of heat advection away from the dayside, so that 
atmospheres experiencing stronger magnetic drag may exhibit reduced infrared phase offsets 
(see, e.g., Fig. 3 of Rauscher & Menou 2011) 

A simple measure of the degree of heat advection and infrared phase offset is provided 
by the ratio r adv / r rad shown in Fig. [2J As T eq increases, r rad decreases steeply (see Eq. [5]) but 
ultimately it is the exponential decline experienced by the zonal wind speed in the strong 
drag limit that leads to the large values of r a d v /T"rad at high temperatures in Fig. |2j One 
would thus expect smaller infrared phase offsets for giant exoplanets with higher radiative 
equilibrium temperatures. It is difficult to evaluate the magnitude of such offsets with the 
simple ratio r a d v /Trad, as opposed to detailed atmospheric circulation models, but we note 
that HD189733b, with T cq ~ 1200 K, still has a measurable offset (Knutson et al. 2007, 
2009a) for a value of r a d v /T r ad ~ 1 (according to Fig. [2] see also discussion in Agol et al. 
2010). 

Interestingly, our scaling law results indicate that the radiative equilibrium temperature 
at which ohmic dissipation peaks is comparable to that at which magnetic drag reduces the 
zonal wind speed by about a half (relative to the case without magnetic drag). A comparison 
between Figures [TJ and H] reveals that this trend holds independently of the planetary magnetic 
field strength. It suggests that, across a range of T eq values, from the start of the magnetized 
regime until the ohmic peak is reached (e.g., ~ 1300-1600 K for _Bdi P = 10 G), zonal wind 
speeds will systematically drop as ohmic dissipation increases. Observationally, this might 
translate into an anticorrelation between the magnitude of infrared phase offsets and the 
degree of radius inflation of hot giant exoplanets. At higher temperatures, beyond the ohmic 
peak, one would expect both zonal wind speeds and ohmic dissipation to drop with T cq . The 
large values of the ratio r a d v /Trad 1) expected at such high temperatures may prevent any 
observable infrared phase offset, but the decline in ohmic dissipation may result in a gradual 
disappearance of the radius inflation phenomenon with T eq , as already mentioned above in 
connection to Fig. 2 of Laughlin et al. (2011). Cowan & Agol (2011) have also noticed that 
hot Jupiters receiving the largest irradiation fluxes (corresponding to T eq > 1900 K) appear 
to be inefficient at redistributing heat to their night sides, a trend which may be consistent 
with the emergence of a strong drag regime. Given the idealized nature of our scaling laws, it 
will be interesting to put these various arguments to the observational test and to interpret 
the results with more detailed, possibly planet-specific atmospheric models. 



-22 - 



6. Conclusion 

We developed scaling laws addressing the magnitude of magnetic drag and ohmic dis- 
sipation in the atmospheres of hot giant exoplanets, as a function of the planetary radiative 
equilibrium temperature. These scaling laws shed light on the nature of the transition from 
an unmagnetized to a magnetized regime of atmospheric behavior, at T cq > 1000 K, and on 
the emergence of a strong drag regime at even higher temperatures. Assuming that poloidal 
field induction remains weak relative to toroidal induction, an issue which deserves further 
scrutiny, our scaling laws suggest that ohmic dissipation eventually vanishes in the limit of 
very strong drag. Clarifying the role of magnetic drag and ohmic dissipation for hot giant 
exoplanets is important and we have suggested that their interplay could manifest in the 
form of an anticorrelation between the amount of infrared phase offset and the degree of 
radius inflation for such planets, thus linking the atmospheric and interior properties of hot 
giant exoplanets in an observationally testable way. 

The author thanks Charles Gammie for providing the original ionization balance routine 
used in this work. This work was supported in part by NASA grants JPLCIT 1366188 and 
PATM NNX11AD65G. 
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